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ABSTRACT 

Weak lensing surveys have become a powerful tool for mapping mass distributions and constraining 
the expansion history of our Universe, but continuum surveys must average over a large number of 
galaxies to average down the ellipticity noise due to the unknown ellipticity and orientation of the 
lensed galaxies. This Letter presents a technique for measuring weak lensing with velocity maps 
that avoids the ellipticity noise. By studying the inherent noise characteristics, we argue this new 
technique could rival or exceed the sensitivity of traditional continuum observations for upcoming HI 
radio surveys. 

Subject headings: gravitational lensing — cosmology: observations — methods: data analysis — 
techniques: radial velocities — radio lines: galaxies 



1. INTRODUCTION 

Weak lensing has become one of the most useful 
techniques for measuring foreground mass distributions, 
and more recently has emerged as one of the most 
promising ways o f constraining the expansion history 
of the universe llMjr^I da-Escudel 119911 iKaiserl Il998t 
iDark Energy Task Forcell200filh 

Weak lensing measurements traditionally work by mea- 
suring the increase in ellipticity and tangential align- 
ment of background galaxie s ijMiralda-Escudel 119911 
iBartehnann &: Schneiderll2001|) . Because the intrinsic el- 
lipticity and alignment of each of the background galaxies 
is unknown, a large number of these galaxies must be av- 
eraged over. This effect is commonly called "ellipticity 
noise," and requires weak lensing measurements to have 
a very high density of background galaxies. 

This Letter presents a technique for measuring weak 
lensing with galactic velocity maps. The use of the pro- 
jected velocity distribution allows weak lensing measure- 
ments based on a single galaxy, and avoids the ellipticity 
noise inherent in broad band lensing surveys. The same 
under lying effect has been previously identified bv iBlair] 
(2002). In this Letter we clarify the underlying physical 
process, how it can be observed, and identify the error 
characteristics. The attributes of the observational er- 
rors are then used to argue that this technique could be 
significant for upcoming HI observations. 

Section[21introduces how lensing can be identified in ve- 
locity maps, and the inherent error characteristics. Sec- 
tion then analyzes the sensitivities, and argues that 
the velocity map technique may be of comparable sensi- 
tivity to broad band observations for upcoming HI radio 
surveys. Whether broad band or velocity map techniques 
are more sensitive is likely to depend on the details of the 
observation and instrument, and the currently unknown 
systematics. 

2. LENSED GALACTIC VELOCITY MAPS 

In standard week lensing measurements the distribu- 
tion of light of the observed galaxy is symmetric around 
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the axis of rotation (in practice and by assumption). 
However, this is not true in velocity maps. The pro- 
jection of the velocity provides a well-defined azimuthal 
structure in the galactic image. It is the azimuthal struc- 
ture in the velocity map and how it changes under shear 
which allows us to measure the lensing signal and elimi- 
nate ellipticity noise. 

In the absence of lensing, a simple disk galaxy will 
produce a cross-like velocity image as shown in panel a) 
of Figure ^ F° r circular orbits the projected velocity 
is simply "^ m " sin 0, where r is the distance from the 
center of the galaxy, m v is the projected distance along 
the major axis, and <f> is the angle between the line of 
sight and the galaxy's rotation axis. This produces a 
maximum in velocity along the m v axis and a zero along 
the mo axis, which are necessarily separated by 90° in the 
unlcnsed velocity map. This simple cross-like structure 
is what would be seen in the absence of a lensing source, 
or equivalently, is the structure seen at the plane of the 
lens. 

To determine the lensed velocity map, we choose axes 
si,S2 aligned at 45° degrees to the axes of the initial 
galaxy. Following Miralda- Escudei l)1991lk the image is 
then stretched by 

[l_(«_ 7 )]-i (1) 
along the s± axis, and 

[1-(k + 7)] _1 (2) 
along the S2 axis, where k is the convergence and 7 is the 
shear along si, S2- The resulting image is shown in panel 
b) of Figure^ Since the vertical (si) and horizontal (52) 
components of the unit vectors mo and m v are scaled by 
differing amounts, the zero and maximum velocity axes 
in the lensed image (m , m' v ) are no longer orthogonal. 

The angle between the velocity axes in the lensed image 
is given by 



a = 2 arctan 
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(3) 



Measuring the angle between the velocity axes in the 
lensed image directly probes the weak lensing parame- 
ters. 
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a) b) 

Fig. 1. — This cartoon illustrates the effects of lensing on a velocity map. Panel a) shows a representative velocity map without lensing. 
The axis of rotation is oriented 60° from the line of sight, with the grey value indicating the observed radial velocity. The image is then 
lensed by s\ & S2 to produce the image in panel b). The orientation of the axis along the radial velocity maximum (m v ) and the zero 
velocity axis (mo) are at right angles in the unlensed image, but have a larger angle a in the lensed image. The observed angle a directly 
measures the weak lensing, as described in detail in the text. In the limit of a shear-only field, a directly measures the component of the 
shear indicated by si &: S2, i.e. the shear component aligned with the bisector of a. (For this example a large value of k = 7 = 0.2, 
equivalent to an isothermal mass distribution centered directly below the lensed galaxy, was used for illustration purposes.) 



The image in panel b) clearly shows the increase in 
cllipticity and tangential alignment of the lensed image 
that are used in standard weak lensing measurements. 
Note however, that while the zero and maximum veloc- 
ity axes are aligned with the major and minor axes of 
the unlensed image, this is no longer true for the lensed 
image. 

To explore the properties of this distortion further, 
we note that a only depends on the convergence k to 
second order, and we can safely neglect this small con- 
tribution for the remainder of this discussion. To first 
order, a then only depends on the component of the 
shear orthogonal to the initial velocity axes mo,m v , or 
equivalently the component of the shear aligned with the 
bisector of a and shown by si,S2 in Figure ^ (Note 
that shear axes are orthogonal under 45° rotations, see 
iCabella fc Kamionkowskl l|2005) for a nice discussion of 
the mathematical properties.) shear along the orthogo- 
nal mo,m v orientation will change the ellipticity of the 
observed image, but will not change the orientation of the 
velocity axes. To first order a measures the component 
of the shear field aligned with the bisector of a. 

The orientation of the shear component measured by 
a is determined by the observed orientation of the lensed 
galaxy. To measure the total shear in a small region thus 
necessitates determining 7 for two nearby galaxies with 
different observed orientations 61,82- If we indicate the 
true shear field as e+,e x (often chosen so e x = 0), then 
the component of the shear measured with each galaxy 
is 

71 = e+ cos 201 +e x sin20i, (4) 

72 = e+ cos 20 2 + e x sin26> 2 . (5) 
These can be easily solved to find 

_ 7 2 sin(26>i) -71 sin(20 2 ) ,„s 
6+ ~ sin(20i -26 2 ) ' U 
If the orientation of the shear is already given by a mass 
model, then only one galaxy is needed to obtain a mea- 
surement. 



An equivalent way of viewing the lensing transforma- 
tion is that the lens sees a flat projected image of the 
galaxy with the characteristic cross-shape in the velocity 
map (mo _L m v as shown in the left-hand panel). The 
observer then sees that image rotated around the S2 axis 
(out of the image plane) by the lens. As the rotation 
around S2 increases, the angle a also increases. Look- 
ing at the lens transformation as a rotation of the initial 
velocity map out of the image plane naturally leads to 
efficient algorithms for detecting the weak lensing signa- 
ture in the observed velocity map. 

2.1. Error characteristics 

The observational applicability of measuring weak 
lensing with velocity maps depends on the error charac- 
teristics inherent to the technique. The three primary er- 
ror characteristics of velocity map lensing measurements 
are: 

Reduced number counts. The most obvious change is 
the reduced number counts. Because the ellipticity noise 
has been eliminated, in principal only two galaxies are 
needed to measure the shear as opposed to hundreds. 
This advantage is offset by the much smaller flux from 
emission lines compared to the integrated luminosity. 
The observational implications come down to the relative 
strength of the line emission, and the galaxy luminosity 
function, and is discussed in more depth in Section 

Uncertainty depends on brightness of the image. More 
importantly, the precision of the lensing measurement 
increases with the brightness of the galactic velocity im- 
age until limited by systematic errors. In standard weak 
lensing observations, bright galaxies are no more useful 
than faint galaxies due to the unknown intrinsic elliptic- 
ity, and all galaxies are weighted the same. With velocity 
maps, the uncertainty in the shear keeps decreasing with 
the uncertainty in the image. This makes bright galaxies 
much more valuable, and changes the way one tries to 
perform the observations (a few bright galaxies may be 
sufficient). 

Reduced systematics. Systematic errors can be more 
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important than raw sensitivity, and velocity map lensing 
observations partially avoid three of the common sys- 
tematics affecting continuum shear measurements. First, 
the velocity distribution is expected to be more circular 
than the distribution of light, and velocity lensing mea- 
surements are unaffected by bright knots in the galactic 
image. Second, the velocity map technique is insensitive 
to intrinsic alignments of the lcnscd galaxies. And finally, 
small errors in the telescope PSF become a second order 
effect instead of first order, and should be slightly less 
problematic for the velocity map technique. 

These error characteristics are encouraging, and in the 
next section we discuss the feasibility of measuring weak 
lensing with velocity maps. 

3. OBSERVATIONAL FEASIBILITY 

Because emission lines are much fainter than the inte- 
grated luminosity, any survey targeting galactic velocity 
maps will necessarily have a lower density of background 
galaxies than an equivalent broadband survey. However, 
as we showed in Section [21 the velocity map survey re- 
quires fewer galaxies due to the absence of ellipticity 
noise. In the end, which survey is more sensitive depends 
on the line-to-continuum luminosity ratio, the luminosity 
function (dN/dS) of the lensed galaxies, and the ability 
of the velocity method to identify the shear signal in low 
signal-to-noise images. 

Radio HI observations are the most promising location 
in the electro-magnetic spectrum for producing sensitive 
galaxy surveys with resolved velocity maps. The HI line 
is bright, the luminosity characteristics are relatively well 
understood, and several upcoming instruments will start 
taking much more sensitive HI observations over the next 
few years. 

The following sections review the relevant observa- 
tional considerations and the sensitivity an HI lensing 
survey might attain. 

3.1. Luminosity functions and SKA galaxy counts 

lAbdalla fc Rawlingsl II2005I) and IBlake et all (l2CKT^) 
have performed extensive studies of the continuum and 
HI radio luminosity functions, and the resulting impli- 
cations for radio cosmology with the SKA. Much like 
the optical galaxy luminosity function, the HI luminos- 
ity function is a power-law (-1.3) with a knee at Afjjj. 
For surveys above the knee, increased sensitivity quickly 
leads to higher number counts. Past the knee, however, 
enhancing the se nsitivity only s lowly increases the num- 
ber counts. (See iZwaan et al.1 l)2003|) for measurements 
o f the HI luminosity function at low redshift.) 

Bla ke et al.l (|2004j) have explored these consequences 
for proposed SKA galaxy surveys. For our purposes 
the most important result is that the number density of 
galaxies is > 10 arcmin -2 in the HI survey and ~ 500 
arcmin for the continuum survey (both assuming 4 
hour SKA integrations). For the HI survey each galaxy 
is detected at > 10er to avoid false identification, thus 
each HI galaxy has sufficient signal/noise to apply the 
velocity technique detailed here. 

3.2. Velocity map sensitivity 

A full simulation of the sensitivity of the velocity map 
technique, incorporating all the instrumental and sys- 
tematic effects that will be important, is well beyond the 



scope of this paper. However, we can construct a illus- 
trative toy model. 

In our model we use a 10 x 10 pixel map, with a mea- 
surement of the HI velocity in each pixel. The velocity 
uncertainty is uniform across the map and no intensity 
information is used (no isophotes, which should further 
constrain parameters). For simplicity a single velocity 
vs. radius template is used (based on the Milky Way), 
and we allow the amplitude of the velocity and the radial 
scale to both vary. Our toy model then has 5 parame- 
ters: the inclination and orientation of the galaxy spin 
axis (both prior to lensing), the amplitude and radial 
scale of the velocity, and the shear. Several hundred re- 
alizations are then created, using the significance of the 
galaxy detection to determine the velocity uncertainty in 
each pixel. 

All five of the parameters are well fit, and normally 
distributed around the expected values. For a single 10cr 
galaxy image, the standard deviation of the shear is 0.02. 
As expected, the uncertainty in the shear decreases lin- 
early with the brightness of the galaxy, decreasing to 
0.007 for a 30<r galaxy, and 0.002 for a IOOct detection. 
Using the galaxy counts and luminosity function of the 
proposed SKA HI survey, each square arcminute would 
have of order 10 galaxies, se veral of which would be sig - 
nificantly brighter than 10cr ( Abdalla & Rawlings 2005). 

These numbers can be compared to the shear uncer- 
tainty of ~ 0.015 per arcmin -2 in the proposed SKA 
continuum survey after using 500 galaxy images to beat 
down the ellipticity noise. While the galaxy counts for 
the velocity technique are much lower, this toy model 
suggests that the shear measurements could equal or ex- 
ceed the continuum observations. 

However, we do not want to over interpret the sen- 
sitivity estimates from this toy model. Realistic values 
for the baseline distribution (PSF), degeneracies due to 
fitting velocity profiles, and systematics related to non- 
circular velocity distributions and warped galactic disks 
all need to be incorporated, preferably based on real HI 
observations. Whether HI or continuum lensing surveys 
are more sensitive will likely depend on these details of 
the observational design and systematics. 

In addition, different applications may favor different 
techniques. For mapping cluster mass distributions, very 
deep HI observations could map to much smaller mass 
scales due to the lack of ellipticity noise. While cos- 
mology surveys are simply concerned with the relative 
signal-to-noise over large volumes. 

It should also be noted that most of the next generation 
radio telescopes, such as the xNTD and SKA, will use 
FX correlators. One feature of FX correlators is that the 
frequency resolution is independent of the bandwidth, so 
the HI and continuum surveys can be performed simul- 
taneously with no sensitivity penalty. 

4. CONCLUSION 

In this paper we have introduced a technique for mea- 
suring weak lensing with galaxy velocity maps. The pri- 
mary advantage of the technique is that it avoids the 
ellipticity noise inherent in continuum weak lensing sur- 
veys. We have also argued that due to the increased 
utility of bright galaxies, the velocity lensing technique 
with planned HI surveys may be competitive with pro- 
posed continuum lensing surveys. 
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Much additional work is needed to evaluate the real- 
world utility of the velocity mapping technique. In the 
future we hope to use VLA and EVLA observations to 
study the systematic limits in real galaxy observations. 
Combined with careful studies of the HI luminosity func- 
tion and instrumental effects of the xNTD and SKA de- 
signs, we hope determine the applicability of this tech- 



nique to future weak lensing surveys. 
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